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Precise measurements are performed on spectral lineshapes of spontaneous Rayleigh-Brillouin scat¬ 
tering in mixtures of the noble gases Ar and Kr, with He. Admixture of a light He atomic fraction 
results in marked changes of the spectra, although in all experiments He is merely a spectator atom: 
it affects the relaxation of density fluctuations of the heavy constituent, but its contribution to the 
scattered light intensity is negligibly small. The results are compared to a theory for the spectral 
lineshape without adjustable parameters, yielding excellent agreement for the case of binary mono- 
atomic gases, signifying a step towards modeling and understanding of light scattering in more 
complex molecular media. 

PACS numbers: 42.68.Ca, 42.65.Es, 42.68.Wt, 51.20.-|-d, 51.40.-|-p 


The spectrum of light scattered in a gas is deter¬ 
mined by the fluctuations of its refractive index [1], or, 
equivalently, by the motion of its molecules. When the 
mean free path between collisions is much larger than 
the wavelength, the scattering spectral lineshape is a 
pure Gaussian, to be understood as a Doppler effect. 
At higher pressures collisional excitations and acoustic 
modes come into play, as was recognized independently 
by Brillouin [2] and Mandelstam [3]. In first approxima¬ 
tion, redshifted and blue-shifted frequency components 
are added to the scattering spectrum with characteristic 
shifts Ah' = Vs k/2'K, with Vg the speed of sound and k 
the size of the scattering wavevector, k/2TT = 2 sin(0/2) /A 
with 9 the scattering angle and A the wavelength of the 
incident light. 

Rayleigh-Brillouin scattering in dilute gases offers a 
sensitive probe of gas kinetics. Understanding the scat¬ 
tered light spectrum involves the linearized Boltzmann 
equation [4], and throughout the years intricate approxi¬ 
mations to the collision integral have resulted in various 
kinetic models for the scattered light spectrum. These 
models may be viewed as a success of statistical physics. 
Still, discrepancies with experiments exist, and the ki¬ 
netic models are generally restricted to simple gases. In 
contrast, the Earth’s atmosphere consists of a mixture of 
gases, each of which explores internal molecular degrees 
of freedom. An important practical application of un¬ 
derstanding such mixtures of gases is in its connection to 
laser light scattering (LIDAR) of the atmosphere [5, 6], 
in particular the ADM-Aeolus mission of the European 
Space Agency for measuring the global wind profile [7]. 

The Tenti model is a well-known theory for the spec¬ 
tral lineshape of scattered light in monomolecular gases 
[8, 9]. The spectrum is determined by the communi¬ 
cation between kinetic and internal degrees of freedom, 
which is characterized by a transport coefficient, the bulk 


viscosity ijf,. The bulk viscosity is a dynamic quantity, 
which is not well known at the GHz frequencies of in¬ 
terest in light scattering. Therefore, 77b was used as an 
adjustable parameter to describe spectral profiles in light 
scattering in both coherent [10-13] and spontaneous [14- 
16] Rayleigh-Brillouin scattering experiments. However, 
the Tenti model is not designed to describe light scat¬ 
tering in mixtures. Nevertheless, applying it to air and 
assuming that air is a fictitious gas with effective values 
for its transport coefficients and molecular mass, yields 
fair agreement with experiments [17, 18]. When devising 
a proper theory for air, one faces the formidable task of 
including both kinetic and internal degrees of freedom for 
several species. 

As a first step in understanding light scattering in more 
complex gases, we will concentrate on mixtures of no¬ 
ble gases using a new experimental setup which provides 
spectra with unprecedented statistical accuracy [19]. We 
will compare these spectra to models with no adjustable 
parameters; the only parameter needed is the atomic di¬ 
ameter, which follows from the well-known value of the 
shear viscosity of the pure noble gas [20]. All our ex¬ 
periments are in the kinetic regime, where the mean-free 
path between collisions is comparable to the scattered 
light wavelength. Interestingly, mixtures of gases with 
very different mass behave in a similar fashion as a gas 
of molecules with internal degrees of freedom. While the 
two components of the mixture briefly can have different 
temperatures, a molecular gas can have different temper¬ 
atures associated with translational and internal degrees 
of freedom. It is the relaxation of these temperature dif¬ 
ferences that determines the scattered line shape. 

The components of our He-Ar and He-Kr mixtures 
have a large mass disparity {Mue/MAr = 0.1002, and 
Mue/MKr = 0.0478. With these different size atoms, 
it is only the heavy ones that contribute to the scat- 
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FIG. 1. In He-Ar mixtures, the light He atom acts as a specta¬ 
tor, it participates in the gas kinetics, but does not contribute 
to the scattered light intensity. The spectra of a mixture of 2 
bar Ar and 2 bar He (a) is very different from a spectrum of 
2 bar Ar (b). The gray line indicates the difference between 
the spectra. 


tered light intensity. The intensity is proportional to the 
square of the optical polarizabilities a, with the ratio 
(ane/cTAr)^ = 1.56x 10“^, and (oHe/cKKr)^ = 5.96x 10“^. 
Thus, the light atoms are spectators, and influence the 
spectral line shape only indirectly through collisions. 
Nevertheless, as Fig. 1 illustrates, their influence can be 
large: adding light atoms to a gas of heavy ones signifi¬ 
cantly changes the shape of the scattered light spectrum. 

A schematic view of the setup for the measurement 
of spontaneous Rayleigh-Brillouin scattering is shown in 
Fig. 2. The light from a narrowband continuous-wave 
laser is scattered off a gas contained in a temperature- 
controlled gas cell. The laser is a frequency-doubled 
Ti:Sa laser delivering light at 403 nm, 2 MHz bandwidth 
and 400 mW of output power. The long-term frequency 
drift was measured with a wavelength meter to be smaller 
than 10 MHz per hour. The scattered light is collected 
at an angle of 90° from an auxiliary focus inside the en¬ 
hancement cavity, in which a scattering-cell is mounted. 
The cell is sealed with Brewster windows. The enhance¬ 
ment cavity amplifies the circulating power delivering a 
scattering intensity of 4 Watt in the interaction region 
[19]. The light that passes through the FPI is detected 
using a photo multiplier tube (PMT) which is operated 
in the photon-counting mode and read out by computer. 
All measurements are performed at room temperature, 
297 ±1 K. 

The scattering angle is determined to be 90 ± 0.9° by 
means of the reference laser beam and geometrical rela¬ 
tions using sets of diaphragms and pinholes present in 
the optical setup. The scattered light is filtered by a di¬ 
aphragm which covers an opening angle of 2°, collected 
by a set of lenses, further filtered by an extra pinhole 
{d = 50 fim) and then directed into a hemispherical scan¬ 
ning Fabry-Perot interferometer, which is used to resolve 
the frequency spectrum of the scattered light. To scan 
the FPI plate distance, the spherical mirror is mounted 
on a piezo-electrical translator, which is controlled by 
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FIG. 2. (Golor online) Schematic diagram of the experimen¬ 
tal setup for spontaneous Rayleigh-Brillouin scattering. The 
laser beam (blue line) is amplified in an enhancement cavity to 
increase the scattering intensity. Scattered light at an angle of 
90° is collimated and directed onto a piezo-scannable Fabry- 
Perot interferometer for spectral analyses, and detected on a 
photomultiplier tube (PMT). 

computer. 

The spectral response S{i') of the Fabry-Perot spec¬ 
trometer was measured in a separate experiment, 
and could be parametrized very well by the formula 
S{i/) = [l-I-4(fz/iZw)^sinc^(7ri//i/FSR)] , with sinc(a:) = 
sin(a;)/x, and where ;zfsr is the free spectral range of 
the etalon, zzfsr = 7553 MHz, and z/w = 139 MHz is 
the Airy-width of the transmission peak. All computed 
model spectra were convolved with 5'(zz), and since the 
free spectral range is relatively small, it is important to 
allow for the periodic nature of S^v). 

The light scattering experiments do not provide an 
absolute intensity, therefore the experimental and com¬ 
puted spectra were normalized such that 7(z/) die = 
1 , where the integral extends over one free spectral range 
(FSR), Vb = i^fsr/ 2. Assuming Poissonian statistics of 
registered photon counts, an estimate of the statistical 
error a{i'i) of measured spectra was obtained from the 
square root of the accumulated photon count Ni at each 
discrete frequency let. It was verified that the fluctua- 
tions iVj' at each Vi were independent. The normalized 
error is then = ]\[-t - de(^'i)]^/o-^(f^i)- If 

the computed line shape model Im would fit the measure¬ 
ment perfectly, then only statistical errors remain and the 
minimum of is unity. The difference between theory 
and experiment will be expressed by y^. 

In the past decades, many ingenious efforts have been 
undertaken to arrive at approximate solutions of the 
Boltzmann equation which are relevant for light scat¬ 
tering. Light scattering involves density fluctuations, 
with the spectrum of scattered light equalling the Fourier 
transform of the density-density correlation function. 
Van Leeuwen and Yip showed that this correlation func¬ 
tion follows from the first moment of the solution of the 
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FIG. 3. Comparison of Ar at p = 2 bar and Kr at p = 1.3 bar 
spectra to the A—moment model for N — 5, 13, 20 and 35. 
The spectra are shown for N = 35, while the convergence 
with increasing number of moments N is demonstrated in the 
top panels. 


linearized Boltzmann equation [4]. 

One such effort is based on the Bhatnagar-Gross-Krook 
(BGK) model, which takes a simple relaxation form for 
the collision integral [21], 

^ + ic-V)f = -aif-fr), ( 1 ) 

with c the molecular velocity, / the position-velocity dis¬ 
tribution function, and fr a reference distribution func¬ 
tion. The latter is determined from the requirement that 
N of its moments are the same as those of the complete 
collision integral for monatomic particles with a re¬ 
pulsive interaction potential [22]. Through increasing 
N, increasingly accurate predictions of light scattering 
spectra can be computed. We show this convergence for 
light scattering in pure samples of Ar and Kr in Fig. 3. 
These experiments are in the kinetic regime, with uni¬ 
formity parameter y of order one. The uniformity pa¬ 
rameter j/ of a simple gas is the ratio of the scattering 
wavelength to the mean free path I between collisions, 
y = l/{kl). The only information that the theory fur¬ 
ther needs is the hard-sphere diameter a of the noble gas 
atom, OAr = 3.66 x 10“^°m and oxr = 4.20 x 10“^°m. In 
the case of Ar, the uniformity parameter is j/Ar = 1.14, 
while for Kr j/Kr = 0.96. The N = 35 model can hardly 
be distinguished from the experimental data. However, 
since the minimum is still larger than 1, a slight but 
signihcant difference between model and experiment re¬ 
mains, although on a relative scale it is < 1%. In fact, 
the agreement is so good that the model might be viewed 
as a benchmark testing for experiments. It should be re¬ 
alized, however, that a monatomic ideal noble gas is the 
simplest system thinkable as there are no internal molec¬ 
ular degrees of freedom. 


The BGK-moment approach was used to develop a 
theory of light scattering in binary mixtures of noble 
gases [23]. There, the reference distribution function fr 
in Eq. (1) was chosen such as to satisfy the principal con¬ 
servation properties of the full Boltzmann collision op¬ 
erators, whilst requiring complete correspondence with 
two-fluid hydrodynamics, i.e. the generalized equations 
of Navier-Stokes and Fourier. Therefore this model pro¬ 
vides a precise transition between the hydrodynamic and 
kinetic regimes [24]. In this sense, the BGK model with 
a judicious selection of the reference distribution fr al¬ 
lows one to include the relevant physical phenomena in 
the kinetic description. This is a great advantage for 
the design of models, as the computation of these vari¬ 
ous contributions is cumbersome. It is expected that a 
model that works well for mixtures, no longer works if 
the density of one of the constituents vanishes, and the 
gas becomes monatomic. This is because in the design of 
the mixture model the focus is on inter-species relaxation 
of temperatures and velocities, and not on the relaxation 
of high-order gradients that determine the shape of the 
spectrum of light scattered from monatomic gases. 

Measured and computed spectra for light scattering in 
He-Ar and He-Kr mixtures are shown in Fig. 4. They 
are characterized by the (partial) pressures and the (par¬ 
tial) uniformity parameters. The definition of the par¬ 
tial uniformity parameters yi = l/(kli) of a binary mix¬ 
ture involves the partial mean free paths li of hard-sphere 
atoms. 


k = 


i=i 


MjMj 


( 2 ) 


with Tii and Mi the number density and atomic mass 
of constituent i, respectively, and = {at + aj)/2 the 
distance between the centers of two spherical particles 
with diameters and aj at the instant of collision. With 
all uniformity parameters of order one, the experiments 
are in the kinetic regime. The computation of the theory 
[23] now also needs the hard-sphere diameter of He, oue = 
2.16 X 10 m, and the atomic polarizabilities of the 
noble gases one = 0.227 x 10“^° Cm^ V“^, aAr = 1-82 x 
10-4° Gm2 V-4, axr = 2.94 x lO’^o Gm^ Vh 

The measured mixture spectra are reproduced well by 
the theory. Although the mixture model is designed to 
represent the relevant interspecies relaxation processes, 
which become less important for asymmetric mixtures, 
the agreement with the measured spectra at 1 bar He 
and 3 bars Ar, and the reverse case, is still excellent. 

Early experiments on mixtures of He and Xe atoms 
were done by Glark [25], but were interpreted in terms 
of hydrodynamics, as a complete kinetic theory was still 
lacking. Light scattering on He-Xe mixtures in a range of 
pressures comparable to ours was studied by Letamendia 
et al. [26], and sizable differences with a kinetic mixture 
model [27] were found. 
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FIG. 4. (Color online) Comparing scattered light spectra of 
mixtures to the mixture model of [23] . (a) Equimolar mixture 
He p — 1 bar, Ar p = 1 bar, uniformity parameters pAr = 
1.63, pHe = 0.54, (b) equimolar mixture He p = 2 bar, Ar 
p — 2 bar, j/Ar = 3.28, yue = 1.08, (c) asymmetric mixture 
He p = 1 bar, Ar p = 3 bar, pAr = 2.96, pue = 1.16, (d) 
asymmetric mixture He p = 3 bar, Ar p = 1 bar, pAr = 
3.59, pHe = 1.00, (e) equimolar mixture He p = 1.3 bar, Kr 
p = 1.3 bar, pKr = 3.28, pue = 0.77, (f) equimolar mixture He 
p = 2 bar, Kr p = 2 bar, pxr = 5.01, pue = 1.17. The lower 
line is the difference between experiment and model. Apart 
from the normalization of the spectra there are no adjustable 
parameters. In (d), the dashed thick (gray) line shows the 
Lorentzian spectral lineshape Id Eq. (3) with Dne-Ar the 
hard-sphere value, Dne-Ar = 1.52 X 10“® m^s”^ The full 
thick (gray) line (partly obscured) shows the purely Gaussian 
spectrum of Ar. 


A striking observation is the narrowing of the mixture 
spectrum in Fig. 4(a,d) when the number of helium atoms 
is increased while keeping the number of Ar atoms the 
same. Collisional narrowing takes place when for Ar, col¬ 
lisions with He atoms become far more numerous than 
collisions with another Ar atom. It stands out in the 
strongly asymmetric mixture of Fig. 4(d). In the kinetic 
case (pAr = C’(l))) the spectral lineshape is a Gaussian, 
/(/) oc exp(—(27r^/A:z;o)^)j with vq the Ar thermal ve¬ 
locity, while in the hydrodynamic limit (pAr ^ 10), the 
spectral lineshape is a Lorentzian: the Fourier transform 


of the fundamental solution of the diffusion equation. 


IdH 


2DHe-Ark^ 

(27r^)2 + (DHe-Arfc2)2’ 


(3) 


with Due-Ar the mutual He-Ar diffusion coefficient. This 
collisional narrowing is known as Dicke narrowing [28]; 
in the hydrodynamic limit it can be explained by the 
decoherence collisions which an atomic scatterer experi¬ 
ences when it straggles diffusively through a dense gas 
of spectators. Remarkably, spectral narrowing is mani¬ 
fested here in non-resonant light scattering. In our case 
the partial uniformity parameter evolves from pAr = 1.6 
at PHe = 1 bar to pAr = 3.6 at pHe = 3 bar. Figure 4(d) 
illustrates that our experimental conditions are still far 
from the hydrodynamic limit, as the experiment is closer 
to the Gaussian than to the Lorentzian lineshape. By de¬ 
sign, our theory for the lineshape embodies both extreme 
cases, and it provides a near perfect reproduction of the 
experimental spectrum. 

We have studied Rayleigh-Brillouin scattering in mix¬ 
tures of noble gases, whose constituents have a very dif¬ 
ferent mass. In all cases, the addition of the light He 
atomic gas has a large influence on the spectral line 
shapes, although He atoms hardly contribute to the scat¬ 
tered light intensity. Density fluctuations in mixtures are 
dominated by the relaxation of temperature and velocity 
differences between the constituent gases. These relax¬ 
ations appear to be captured adequately by the mixture 
model of [23]. This is the first time that these predictions 
are tested in the kinetic regime. The model contains no 
adjustable parameters, and reproduces all experiments 
excellently. The present successful comparison of model 
and experiment marks a step towards a description of 
light scattering in mixtures of molecular gases (such as 
air), and emphasizes the need to account for interspecies 
relaxation of temperatures of both translational and ki¬ 
netic degrees of freedom. 
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